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Abstract 


Simplified models, performance correlations and rigorous simulation models are identified 
as alternative techniques for the prediction and characterisation of the performance of natural- 
circulation solar energy water heaters. The features of particular examples of each approach 
are described. In their different appropriate contexts, each technique is shown to provide good 
agreement with measured system behaviour. A more general role in practical system design 
is foreseen for rigorous detailed simulation. Models have been used primarily for generic 
research optimisations. © 2000 Elsevier Science Ltd. All rights reserved. 


1. Introduction 


The heat transfer fluid circulates from the solar-energy collector to the hot-water 
store in a thermosyphon solar-energy water heater by buoyancy forces. Neither a 
pump nor its associated control sensors and actuators are required to circulate the 
heat-transfer fluid. Such solar water heaters thus have both lower initial and running 
costs compared with equivalent forced-circulation units. Furthermore, thermosyphon 
solar-energy water heaters can perform as effectively as forced-circulation units [4]. 

The performance of natural-circulation solar-energy water heaters may be pre- 
dicted analytically by three broad alternative approaches: 
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Nomenclature 

A area (m°) 

C mean specific heat capacity (J kg ! K~') 

D diameter of the collector flow channel (m) 

F collector's efficiency factor 

f solar fraction 

H total daily insolation upon collector (Jm ?) 

h height (m) 

I global insolation falling upon the plane of collector (Wm ?) 
L length of flow channel (m) 

M mass of water (kg) 

m mean daily thermosyphonic mass flow rate (kg s) 

N number of times the water contained in the store is circulated 


through the collector during the insolation period 


n integer 

QO mean rate of delivery of useful heat from the collector during the 
insolation period (W) 

Q total heat delivery (J) 

s collector-plane's inclination to the horizontal (Degrees) 

T temperature (K) 

t time (s) 

U steady-state heat loss co-efficient (Wm ? K ') 

X,Y Z dimensionless groups defined by Eqs. (26), (27) and (28), 
respectively 

a absorptance of the absorber surface 

B coefficient defined by Eq. (32) 

AT temperature difference between the fluid at the inlet and outlet of 
the collector (K) 

1] efficiency 

0 duration of insolation period (s) 

T transmittance of collector glazing 

6 L-T 

Subscripts 

a of the ambient environment 

c for the collector 

d refers to the downcomer pipe 

e effective 

f for the heat-transfer fluid 

1 refers to collector inlet 

I refers to water withdrawn from store 
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mean value 

refers to number of passes 

for the collector plate 

heat-removal factor 

refers to the system as a whole 

for the hot-water store 

refers to the upriser pipe 

for the volume of water to be heated 

refers to the fully-mixed store, immediately prior to the insolation 
period 


=N E NQ% WDB 


Superscripts 


f denotes absorber plate efficiency factor 


- denotes average value of the variable 


e simplified models, providing average or cumulative diurnal predictions 

* correlation of performance characteristics from either simulation or monitoring of 
generic systems, and 

* rigorous dynamic simulation models using detailed numerical analysis, providing 
hour-by-hour or more frequent performance predictions 


Simplified models 


A multiple-pass solar water heater is defined as a system in which the volume of 
water in the store is circulated through the collector many times each day. It has 
been observed in numerous experimental studies that, for a multiple pass thermosy- 
phon solar-energy water heater, the difference in the fluid temperatures at the inlet 
and outlet of the solar-energy collector is constant [5]. Whether a particular system 
is single or multiple pass depends on the system’s geometry and flow resistance as 
well as on the prevailing insolation level. An almost constant temperature differential 
across a collector (concomitant with multiple-pass behaviour) ensues, provided either 
a critical height difference between the collector and the elevated hot-water store for 
a fixed rate of heat input, or a critical insolation level for a fixed system geometry 
has been exceeded [6]. 

The exhibition of an almost constant temperature differential across the collector 
has been used to develop a particular simplified analysis — the ‘AT model’ — for 
the determination of daily system performance [5]. The advantage of this particular 
approach is that there is no necessity to determine the frictional flow losses within 
the fluid circuit nor to perform a time-wise analysis of the temperature variation 
around the thermosyphon loop to ascertain the buoyancy pressure head. 
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To develop any simplified model, many assumptions are required. For the AT 
model’, these are: 


1. the temperature differential, AT, across the collector remains constant throughout 
the insolation period. The duration of periods at the start and end of each insolation 
period during which the system does not exhibit this behaviour are assumed to 
be much shorter than the insolation period. 

2. diurnal mean insolation and ambient temperatures can be employed in the analysis 
without introducing significant errors in the calculated values of the net daily 
heat gain. 

3. at the start of the insolation period, the hot-water store is fully mixed at a tempera- 
ture 7;. During the period of insolation, the mean store temperature rises linearly 
with time to a value T, at sunset. No hot water is withdrawn from the system 
during this period. The latter assumption implies that the system is intended to 
produce hot-water for withdrawal either that evening or early the subsequent day. 

4. all fluid properties are independent of temperature. 


For no mixing in the hot-water store, the water temperature in the store after n passes 
of the full volume of the store through the collector is: 


T,=T,+n AT (1) 


This is also the inlet temperature, 77; for the (n+1) pass. If the total number of passes 
is N, then at the end of the day 


T,=T,+N AT (2) 
"EIL N 

1-5) ren (or (3) 
srt (4) 


Alternatively, when the store is fully mixed, the rate of increase of the store tempera- 
ture is: 


dT m 
aed 6) 
thus 
m 
T(th=T,+AT+ (5) (6) 


At the end of insolation period, 
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m 
T.=T,+AT (s) (7) 


The inlet temperature at the collector is assumed synonymous with the tank tempera- 
ture, throughout the insolation period. The diurnal average fluid inlet temperature at 
the collector can be found by taking the diurnal average of the tank temperature, i.e., 


2 1 m 
T,;=T, «(S)are() (8) 


Thus the mean fluid temperature over the whole day at the inlet to the collector is 
given by Eq. (4), irrespective of the degree of mixing in the store. The useful rate 
of heat provided by a solar energy collector can be described by: 


Q- FaA Iro — UT;;—T,)] (10) 


where Fg, the heat removal factor is given by 


F= ( exp UA) (11) 


mC; 


The rate of energy delivery by the collector is: 
Ò=rc AT (12) 


Substituting from Eqs. (4), (10) and (11) into Eq. (12) gives: 


A. =c; U.AT 


= —log.4 1 
T+T, 
Ita—U 2 “J| I-T, 


m U.F' 
From the energy balance for the whole system, over the full insolation period the 
average collector plate temperature is: 


TE (; Er) 14 
P FR Fr ü Mita 2 ( ) 


for which it is assumed that the heat losses from the up-riser and down-comer are 
negligible, thus: 


(13) 
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T, (nc, - T marc 
T= 


i U 
m,C,— vot 


HT ,U A0 (15) 


The system’s diurnal heat-gain efficiency is defined as 


M;CKT.-T) 


il IOA. (16) 


The critical minimum vertical height of the down-comer which will prevent circu- 
lation occurring at night is given by [5]. 


1 
hum cp Ta Tosin s+ (L.-T) i 


If, for the particular installation, the assumption that, at night, 


TAT, 


T,—-T, Ta=T,+5 (in°C) and T,= 2 


adequately reflect reality, then Eq. (17) becomes 


———— c €— h (18) 
d I : c 2 T 
2 


Using initial estimates for T. and Tp, U, is calculated. Using Eq. (14) a new value 
of T, is calculated, this process is repeated until convergence. A new value for T, 
is found from Eq. (15) and, using the convergence value of T,, an iterative procedure 
is used to obtain a solution for T.. If mutually-incompatible parameters are introduced 
into the procedure, then no solution will be found. This is because Eq. (13) cannot 
be solved as: 


(19) 


N TaT. 
.ar>ta-u4| - 2r. 


When compatible values for the overall collector heat loss co-efficient are chosen, 
the condition expressed in Eq. (19) will not arise; for example if the value of U, is 
large then a small AT would be expected. 

Theoretical predictions from the AT model have been compared with measured 
performance data from two identical thermosyphonic solar-energy water heaters [6]. 
Each 1.4 m? collector consisted of single-glazed aluminium flat plate absorbers, 
painted matt black: their fourteen flow channels were of 25 mm internal diameter. 
They were inclined at 60° to the horizontal and were oriented facing south. The 
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storage tanks were each 1.6 m tall and 0.3 m in diameter, giving an aspect ratio of 
5.3 and capacity of 0.1125 m°. In each system the upriser discharged into the tank 
at a distance of 30 mm below its top, and the downcomer was connected at the level 
of 30 mm above the tank base. The data acquired included the time histories of 
temperatures at various heights in the collector store and the ambient temperature 
using copper/constantan thermojunctions. The global insolation on the collector sur- 
face was measured with a Moll-Gorczynski pyranometer, mounted in the same plane 
as the collectors. During the tests, no water was withdrawn from the store. Compari- 
sons of the predicted and observed values for the mean temperature of the stored 
hot water are shown in Fig. 1. These results, obtained for a wide range of ambient 
and initial system conditions, show good agreement. The model could not be applied 
to the system for those days during the test period with a final mean store temperature 
of <38°C, as the collector temperature differential did not exhibit a constant value 
throughout the majority of the insolation period. 


indicating possible error of +1°C 


52 t, 
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» 
48 t + 
4 
"^. 
of 
46 / 
f- 


i 74 
zd 


[ 1 1 r J 


l = 
40 42 44 46 48 50 52 
Predicted mean temperature of the stored hot water 
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Fig. 1. Comparison of the experimentally-observed final mean temperature of the stored hot water with 
that predicted by the ʻAT’ simplified model. No water withdrawal from store. 
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3. Correlation of performance characteristics 


In this approach, using simplifying assumptions, an analysis of the thermal behav- 
iour is undertaken to identify dimensionless groups. The developed correlations 
between these characterising groups describe long-term system performance. As the 
simplifying assumptions do not apply in all cases, the functional relationship between 
the dimensionless correlating parameters obtained analytically is not that observed 
in reality. The actual correlation between the identified grouped parameters is 
obtained either by comparison with the output of a rigorous numerical simulation 
model or data from long-term system performance monitoring. 

In a particular approach to the correlation of the daily performance of indirect 
thermosyphon solar-energy water heaters [2], a transient heat balance was undertaken 
on a generic solar-energy water heater (see Fig. 2), with the following assumptions: 


Losses to the Water drawn off 
atmosphere Temperature, T, 
(UA), (T, -T, — 
AMBIENT 
ENVIRONMENT Insulated | Net thermal load 
M,C, (T, - To) 
Net heat transfer A 
from collector to 
store Z — 
Fay A, (T &1- U (TT) Make-up water at 


Temperature, Tm 


Collector losses 
to atmosphere Insulated 
Fay AL oTa) downcomer 
Flat-plate Collector 
solar-energy 1 
colector rud gains 
FayAT OL ATTIC SPACE 


Air gap 


=> _ Indicates direction 
of thermosyphon 
flow during heat 
collection 


— —v-- Indicated components 
of energy balance 


Fig. 2. Holistic heat transfer mechanisms for a generic thermosyphon solar energy water heater. 
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(i) the collector and storage tank both have the same mean diurnal temperatures, 
(i1) insolation and ambient temperature are assumed to remain invariant over the 
day at their respective mean measured values, 

(iii) the total draw-off of water from the store takes place as a single event at the 
end of the insolation period, 

(iv) water is drawn off at the mean temperature of the store and replenished at 
the mains water supply temperature, 

(v) tank and pipe losses are negligible, and 

(vi) the store is initially at the mains water supply temperature at the onset of the 
insolation period. 


An instantaneous heat balance on the system gives, with assumption (i) 


M EE LT Tat, 20 
P dt S AL (000.5 ( F wl ( ) 


Using assumptions (ii) and (v), after making the substitution and rearranging, Eq. 
(20) becomes, 


dó 1 FRA UQ | FaA (10).H 


| 21 
dt M,Cp M,C,0 en 
With assumption (vi) as a boundary condition, the solution to Eq. (21) is 
FRA (cO) H-FpA-UO(T,,-T,) FRAU 
= - 1 22 


The difference in the mean store temperature between the beginning and the end of 
the insolation period is therefore 


| | FA (10), H—FgA,UO(T,,-T,) FRA.U,@ 

? zd FAUO IL e M.C, ) s 

From assumptions (iii) and (iv), the solar fraction can be expressed as: 
M, C.($— 

jeg A e (24) 

Substituting Eq. (23) into (24) gives: 
SQ _ FrA(ta).H (25) 
Frys Uo F,A,UO(T,,-T,) 
MCa aTa) 1—e- M,C 


from which three dimensionless groups can be defined as: 
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For the particular thermosyphon solar-energy water heater described, inspection 
of the dimensionless groups X, Y and Z indicates a linear relationship, shown graphi- 
cally in Fig. 3, between the grouped parameters is applicable to A(ta), Fp, Ur, Ac 
and M; which defined the system configuration and thermal characteristics. The per- 
formances of buoyancy-driven systems also depend on other factors such as the 
position of the store relative to the collector. The influence of these other factors 
may be encapsulated in the relationship between the store and collector temperatures, 
that is the validity, or otherwise of assumption (i). An initial investigation indicated 
a constant temperature difference between the daily mean collector (T,) and store 
(T7) temperature over the insolation period. Fig. 4 shows a plot of the daily values 
for T, against Ts for the month of August indicating a linear relationship. The best 
straight line fit through the data gave: 


T,—1.00 Tj,—7 —1.41 (31) 
with a correlation coefficient of 0.97. Rearranging Eq. (31) gives: 
Tr—T,—1.41 


Thus, more generally 
T,—T.—B (32) 


Introducing into the analysis such a constant temperature difference B, between the 
store and collector results in the dimensionless group Y being redefined as, 


600r 


s 
2 


[m] 


4^. Best fit straight line 
through data 


ü 


o 
T 


DP 


X/[1-exp(-Z)] (DIMENSIONLESS) 
8 
o 


Jd. Measured daily values 


-200 + 
a 
0 «8000 4000 
Y/Z (DIMENSIONLESS) 
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Fig. 4. Correlation of the mean collector and hot water store temperatures, each averaged over single 
periods of insolation. 


_FavAd(Tæ)H+U p0] 
MCT nTa) 


Using this modified value of Y results in a correlation with a more universal appli- 
cation potentially. The performance of any system can then be determined from 
several daily measurements of D. Experimentally-determined values of B would 
enable X, Z and a modified Y parameter to be used to predict the daily performance 
of the system under the applied conditions with a universal correlation. Alternatively, 
daily values of B could be generated using a high-level rigorous simulation model 
and correlated against the height of the store relative to that of the collector. The 
latter approach is employed by [7] to develop a single sizing-chart nomogram. 


Y 


(33) 


4. Rigorous dynamic simulation models 


In a rigorous model, the coupled energy and momentum equations which describe 
the transient behaviour of a natural-circulation solar-energy water heater are solved 
numerically. The main features of a representative example of a dynamic simulation 
model [3] are: 


(i) a finite difference, transient heat transfer analysis is applied to the liquid circul- 
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ating through all the components (i.e., collector, upriser, hot-water store and 
downcomer) of the thermosyphon loop, 

(ii) the thermal capacitances of both the collector plate and the cover are included, 
(iii) the density, specific heat, viscosity, conductivity and Prandtl number of the 
circulating liquid are represented as second-order polynomial functions of tem- 
perature, 

(iv) all the heat-transfer coefficients are temperature dependent and based on the 
ambient and mean component temperatures and are updated at each time step in 
the numerical calculation, 

(v) the store model includes buoyancy-induced mixing between stratified layers 
which occurs when warmer fluid is introduced below a cooler layer, 

(vi) a transient momentum equation is employed, 

(vii) friction factors are calculated using correlations appropriate to both the non- 
isothermal thermally-destabilised low Reynolds number flow and isothermal 
developing laminar flow in straight pipes, 

(viii) empirically determined laminar heat-loss coefficients for the pipe bends 
are included, 

(ix) time variations of insolation, ambient temperatures and hot-water withdrawal 
encountered by an installed system during operation [8], were used as inputs to 
the simulation, 

(x) the transmission of the glass collector cover is determined as a function of 
the sun hour angle, and 

(xi) a two-dimensional finite-difference model for the flat-plate solar-energy col- 
lector was employed. The main advantages of this over previously used pro- 
cedures are 

(a) thermal masses and temperature profiles for each of the fluid, collector plate 
and glass cover can be described individually; and 

(b) fewer simplifying assumptions are made about both the heat transfer processes 
within the collector plate and between the collector plate and the fluid. 


The equations are expressed in finite difference form and solved simultaneously using 
the appropriate boundary conditions. In this model an implicit method of solution 
was employed in which N simultaneous equations are set up for the N nodes. and 
equations were solved using a Gauss-Seidel iterative method. For this method, the 
solution is stable unconditionally and the size of the time step is only limited by the 
required accuracy of the prediction of performance. 

Validation was undertaken in an indoor thermosyphonic solar-energy water heater 
test apparatus, the measured global insolation was simulated experimentally by reg- 
ulating, via a computer, the electrical power supplied to the heater. Similarly, the 
thermal load on the system was reproduced experimentally by the controlled opening 
and closing of the solenoid valves giving a draw-off pattern equivalent to that pro- 
duced by the occupants of the house. The measured insolation, draw-off pattern and 
ambient and mains temperatures were used as input data to the performance simul- 
ation. Measured values of the global insolation and water draw-off patterns were 
obtained from an occupied dwelling in the south-east of England to which a thermo- 
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syphon solar-energy water-heater had been installed and monitored [7]. Fig. 5 shows 
the variation of the experimental and predicted mean store temperatures during the 
day. The predicted mean store temperatures agree to within 0.7?C of the correspond- 
ing measured values over the period of time from 00.00 hours to approximately 20. 
10 hours, when a 12 kg water draw-off was made. After this draw-off, the predicted 
mean store temperature fell to 1.3?C below the measured value, i.e. the predicted 
temperature of the water drawn off from the top of the store was higher than the 
actual temperature of the draw-off water. 

Good agreement is shown, as illustrated in Fig. 6, between the measured and 
predicted values for both the collector's inlet and outlet temperatures over the dur- 
ation of the simulated insolation period. Discrepancies between these sets of values, 
particularly at the end of the insolation period, are due predominantly to the uncer- 
tainty in predicting the frictional flow losses at the very low Reynolds number flows 
encountered during reverse circulation. The predicted flow rates shown in Fig. 7 
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Fig. 5. Diurnal variations of experimental and predicted, via a rigorous simulation model, mean store 
temperature of a thermosyphon solar energy water heater. 
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Fig. 6. Experimentally-observed inlet and outlet temperatures of a flat-plate collector in the flow circuit 
of a thermosyphon solar-energy water heater compared with those predicted from a rigorous simulation 


model. 


correspond to within 7% of the measured values over the insolation period. Flow- 
rates were not measured during periods of reverse circulation as no signal was pro- 
vided by the flowmeter below a positive flow of 10 * kg s™'. 

The large fall-off and subsequent recovery in the measured flow rate immediately 
after water had been drawn off from the store was predicted correctly via the math- 
ematical model. This reduced flow rate occurs as the drawn-off water is replenished 
by cooler water (at the mains supply temperature) introduced at the base of the store. 
This cooler fluid enters the collector and the resulting increase in fluid density in 
the risers and upriser pipes reduces the pressure differential between the fluid in 
these components and that of the store. A subsequent fall in the flow rate ensues. 
The reduced flow rate caused the fluid temperatures in the collector to increase; the 
fluid density falls and a higher rate is restored rapidly. The total thermal output due 
to the 7 draw-offs made during the day was measured as 4.27x10 * J. The corre- 
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Fig. 7. Diurnal variation of experimentally-observed thermosyphonic mass flow rates and those predicted 
from a rigorous simulation model. 


sponding predicted value was 4.54x10 ? J. This larger heat output can be attributed 
directly to the higher degree of stratification in the store predicted by this particular 
simulation model. Rigorous dynamic simulation models have been used for both 
direct [2] and indirect [9] thermosyphon solar water heaters, the latter have heat- 
exchangers in the store through which an anti-freeze solution is circulated to the 
collector. Detailed rigorous simulation models have been employed to determine the 
optimal concentration for an anti-freeze solution [8]. A dynamic simulation model 
with a three-dimensional representation of the hot water store has been employed to 
analyse the effect of tank geometry on thermal stratification [1]. 


5. Conclusions 


The three generic approaches discussed for the analysis of thermosyphon solar 
energy water heaters are summarised in Table 1. The simplified approach described, 
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Table 1 
Salient features of alternative generic approaches to the analysis of thermosyphon solar-energy water heat- 
ers 


Roles Limitations Advantages 
Simplified e System sizing * Only apply to the range of systems and * Small input data 
models operating conditions for which the set 


simplifying assumptions are valid. 

* Otten require experimentally-determined 
information only obtainable after the 
construction of the unit. 


Correlation * Long-term * Cannot be applied reliably to those * Small input data 
performance performance systems for whose dimensions and set 
characteristics prediction climatic conditions, a correlation has * Wide-range of 
Design charts not been determined. systems can be 
Economic considered 
optimisation * Designer 
orientated 
Rigorous * Engineering * Require increased computing power * Accurate 
simulation optimisation * Large input data set transient 
models * Establishing simulation 
long-term * All system 
performance configurations 
* Determining can be included 
validity of 
simplified 


models 


i.e. the AT model, is very quickly executed. Its application is limited by the need 
for the mean diurnal insolation intensity to exceed a, usually unknown, threshold 
level for a given system geometry. The principal generic shortcoming of such simpli- 
fied models is that there is no prior way of determining the additional performance 
parameter, in this case the collector temperature differential AT, knowledge of which 
facilitates the simplification of the analysis, without performing either an experi- 
mental measurement or a rigorous simulation. However, for the particular case of 
the AT model, an analysis of the effects of changes in the AT value employed 
revealed that though the fluid mass flow rate is very sensitive to this factor, the final 
fluid temperature is much less affected. 

The performance correlation described enables the long-term monitored perform- 
ance of an indirectly heated thermosyphon solar-energy water-heater to be related, 
on a daily basis, to the thermal demands made on the system and the prevalent 
meteorological conditions. The correlation exhibits a high degree of linearity over a 
wide range of climatic conditions, and draw-off patterns. 

Rigorous dynamic simulation models have required ‘expert’ users. However rapid 
and continuing advances in graphical user interfaces and pre- and post-processes 
together with the ready availability of computing capacity may render the bespoke 
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design of large-scale non-domestic thermosyphon solar water heaters practicable 
using such models as tools. 
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